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We present a quantum network approach to the treatment 
of thermal and quantum fluctuations in measurement devices. 
The measurement is described as a scattering process of input 
fluctuations towards output ones. We present the results ob- 
tained with this method for the treatment of a cold damped 
capacitive accelerometer. 



I. NON IDEAL QUANTUM MEASUREMENTS 

Active systems are fundamental elements in high pre- 
cision measurements. Amplifiers are used cither for am- 
plifying the signal to a macroscopic level or to make the 
system work around its optimal working point with the 
help of feedback loops. With techniques such as cold 
damping, it is possible to manipulate actively the fluc- 
tuations and to reduce the effective noise temperature of 
the devices well below the operating temperature. The 
analysis of sensitivity limits in these devices rises many 
questions related to fundamental processes as well as ex- 
perimental constraints. How far is it possible to reduce 
the measurement temperature? How are these process re- 
lated to the fluctuation dissipation theorem? Are there 
quantum limits to this noise reduction associated with 
Heisenberg inequalities? How do the experimental con- 
straints interplay with the fundamental limitations of the 
sensitivity? 

The aim of the present paper is to address these ques- 
tions with quantum network theory. This approach pro- 
vides a rigorous thermodynamical framework able to 
withstand the constraints of a quantum analysis of the 
measurement. In the same time, it makes possible a re- 
alistic description of real measurement devices. Ther- 
modynamic and quantum fluctuations are treated in the 
same footing. The measurement process is described as 
a scattering process allowing for a modular analysis of 
real quantum systems. Active systems such as the linear 
amplifier or the ideal operational amplifier are described 
in this framework. Here, the approach will be illustrated 
by analyzing the sensitivity of a cold damped capacitive 
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accelerometer developed for fundamental physics appli- 
cations in space . 

We first present the analysis of passive electrical sys- 
tems in term of quantum networks. Then, we use this 
approach to present the quantum analysis of an opera- 
tional amplifier working in the ideal limit of infinite gain, 
infinite input impedance and null output impedance. In 
the last section, we illustrate the theoretical framework 
with the example of a cold damped accelerometer. 



II. COUPLING WITH THE ENVIRONMENT 

Relations between fluctuations and dissipation have 
first been discovered by Einstein who studied the viscous 
damping of mechanical systems Another important 
application was the study of Johnson-Nyquist noise in 
resistive electrical elements IBJ. This classical result was 
extended to take into account the quantum statistical 
properties of fluctuations A general approach of 

these relations was widely studied in the framework of 
linear response theory MS . 



A. Dissipation and Fluctuations 

A first insight into the physical effect of the coupling of 
an electrical circuit to the environment is provided by the 
analysis of an antenna in an electrical resonator. When 
a current flows through the antenna, electromagnetic ra- 
diation is emitted and the resonator energy decreases. 
As far as the electric circuit is concerned, the effect of 
the antenna is the same as a resistance. The antenna 
is also able to detect electromagnetic fields. An incom- 
ing wave puts into motion the electrons in the antenna 
and causes an electrical current to flow in the circuit. 
For thermal radiation, the detection radiation leads to a 
random current which brings the electrical oscillator to 
thermal equilibrium. In the high temperature limit, it 
leads to the usual thermodynamic \ksT per degree of 
freedom, with ks being Boltzmann constant and T the 
radiation temperature. In the zero temperature limit, the 
detected field corresponds to the vacuum fluctuations of 
the electromagnetic field and the induced energy of the 
oscillator is the zero point energy ^hujQ, with oj the res- 
onance frequency of the oscillator. 
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FIG. 1. Representations of a resistance R. (a) Thevenin 
representation with a voltage noise generator U n . (b) Model 
with a semiinfinite line and propagating fields a m and a 

In figure [l] are depicted two representations for a resis- 
tance R. Figure [l]a corresponds to the Thevenin repre- 
sentation with the noise source represented as a voltage 
noise generator U n . The relation between the current I 
and the voltage V is 

U = RI + U n (1) 

Figure |l| b corresponds to a model that originates from 
Nyquist's analysis 0. It consists in a semi infinite coaxial 
line of characteristic impedance R. The solution of the 
propagation equations in the line may be written as the 
sum of two counterpropagating fields 7 out and I ln 

I(x,t) =/ out -7 in (*-~) 

U (x, t) =R (l out (t + ^ ) + / ln (t - ^)) (2) 

At the end of the line, we deduce the following relations: 
U = RI + 2RI in = RI + U n 

jout _ J _j_ Jin (3) 

The first equation corresponds to the relation (Q) and 
leads to the identification of the noise as the input cur- 
rent 7 m . The second equation describes the output fields 
jout em jtt ec i back to the line. This output field may be 
used either to feed other elements of the system or to 
perform a measurement by extracting information on the 
system of interest through a line considered as the detec- 
tion channel. 



B. Treatment with quantum fields 

In an infinite line, current and voltage may be treated 
as quantum fields propagating in a two dimensional 
space-time. Throughout the paper, we will consider that 
a function / is defined in the time domain (notation 
/ (i)) or in the frequency domain (Kubo's notation / [lo]) 
and that these two representations are related trough the 
Fourier transform with the convention of quantum me- 
chanics 

/(*) = / ^/He^ (4) 

The electronics convention may be recovered by substi- 
tuting j to —i. 



Free field operators a m and a out can be defined as the 
Fourier components of J m and I out 



I(x,t) = 



U(x,t) = 




They are normalized so that they obey the standard com- 
mutation relations 

[a in [u] , a in [u/]] = [a out [u] , a out [to'}] 

= 2tt 5 (lo + lo') e (lo) (6) 

where e (lo) denotes the sign of the frequency to. This 
relation just means that the positive and negative fre- 
quency components correspond respectively to the anni- 
hilation a u and creation aj, operators of quantum field 
theory 

a in [u] = a^B (lo) + alj {-to) (7) 

9 (lo) denotes the Heavyside function. 

To characterize the fluctuations of these noncommut- 
ing operators, we use the correlation function defined as 
the average value of the symmetrized product. With sta- 
tionary noise, the correlation function depends only on 
the time difference 

{a in (t)-a in (t')) =*Z(t-t') 

(a in M-a in K]) =27r£(a,W) o*[w] (8) 

The dot symbol denotes a symmetrized product for quan- 
tum operators. 

In the case of a thermal bath, the noise spectrum is 

aa[ J exp^- 1 2 2 2k B T a 1 ' 

One recognizes the black body spectrum or the number 
of bosons per mode for a field at temperature T a and a 
term i corresponding to the quantum fluctuations. The 
energy per mode will be denoted in the following as an 
effective temperature a 

k B e a = n m az h = coth 1\± ( io) 

2 2k B ±a 

In the high temperature limit the classical energy for an 
harmonic field of k B T a per mode is recovered. In the 
low temperature limit, the energy corresponding to 
the ground state of a quantum harmonic oscillator is ob- 
tained. Note that the term i corresponding to the zero 
point quantum fluctuations was added by Planck so that 
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the difference with the classical result ksT a tends to zero 
in the high temperature limit JIo[| . 

These results are easily translated to obtain the ex- 
pression of the Johnson Nyquist noise power 



a UnUn [w] = 2Rh \u\ aZ [u] = 2Rk B Q a 



(11) 



Our symmetric definition of the noise power spectrum 
leads to a factor 2 difference with the electronic conven- 
tion where only positive frequencies are considered. 



C. Quantum networks 

The elementary systems described up to now as well 
as more complex devices to be studied later in this paper 
may be described by using a systematic approach which 
may be termed as "quantum network theory" . Initially 
designed as a quantum extension of the classical theory 
of electrical networks , this theory was mainly devel- 
oped through applications to optical systems |p^| , |l3| . It 
has also been viewed as a generalized quantum exten- 
sion of the linear response theory which is of interest for 
electrical systems as well [Q. It is fruitful for analyz- 
ing non-ideal quantum measurements containing active 
elements 15 lq|. 

In this quantum network approach, the various fluc- 
tuations entering the system, either by dissipative or by 
active elements, are described as input fields in a number 
of lines as depicted on |l| b. 
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FIG. 2. Representation of an electrical circuit as a quan- 
tum network. The central box is a reactive multipole which 
connects noise lines corresponding to the fluctuations entering 
the system, either by dissipative or by active elements. For 
example, the upper left port n with voltage U n and current 
I n is connected to a line of impedance R n with inward and 
outward fields ajf and a„ . 

We first consider a passive linear network built with 
resistances and reactive elements like capacitances or in- 
ductances. Each resistance R n is modeled as a semi- 
infinite coaxial line a n with characteristic impedance R n . 
The voltage U n and current I n associated with the re- 



sistance are the inward and outward fields a" 
evaluated at the end of this line 



and 



R-»y5H ( , 



U = R" 



i n\ 



(a out + a in ) 



(12) 



Here, X = I, U, a ln , a out denotes the column vector with 
components X n and R is the diagonal matrix formed 
with the characteristic impedances R n . 

Input fields corresponding to different lines commute 
with each other. For simplicity, we also consider that 
the fields entering through the various ports are uncorre- 
cted with each other. The interaction with the reactive 
elements is described by a reactive impedance matrix Z 



U 



Z I 

z 



(13) 



The whole network is then associated with a scattering 
S matrix, also called repartition matrix Jl7t , describing 
the transformation from the input fields to the output 
ones 



= S a m 
RiZR 



R iZR 



(14) 



The output fields a out are also free fields which obey the 
same commutation relations [}] as the input ones. In other 
words, S matrix is unitary. In the case of the passive 
network, this property is an immediate consequence of 
the reactive nature of the impedance matrix Z 



R sZtR s_i -R-sZR-'-l 
R iZtR 3+1 -R iZR-5 + 1 



(15) 



More generally, the unitarity of the S matrix is re- 
quired to ensure the quantum consistency of the descrip- 
tion. In the following section, we will make use of this 
property to deduce general properties of amplifiers. 



III. FLUCTUATIONS IN AMPLIFIERS 

Quantum noise associated with linear amplifiers has 
been the subject of numerous works. In the line 
of thought initiated by early works on fluctuation- 
dissipation relations, active systems have been studied in 
the optical domain when maser and laser amplifiers were 
developed fl8|-|20[| . General thermodynamical constraints 
impose the existence of fluctuations for amplification as 
well as dissipation processes. The added noise determines 

and 



the ultimate performance of linear amplifiers 21 22 
plays a key role in the question of optimal information 
transfer in optical communication systems |25| , |2"4j| . 

We first consider the amplification of a field, for ex- 
ample in long distance telecommunication systems with 
repeaters. 

The amplification of the field a m with a phase insensi- 
tive gain G is given by the following equation 



Ga ln + B" 



(16) 
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where B ln is a noise added by the amplification. The 
gain G may be frequency dependent. The commutator 
of the output field is then 



[a° ut M,a° ut M] = \G\ 2 [a h 



+ [B in [<j] , B in [a/]] (17) 

The unitarity of the input output transformation and 
the preservation of the commutation implies a non zero 
commutator for B m 

[B in [lu] , B in [u/]] = [a in [lj] , a in [a;']] - [a out [w] , a out [a/]] 
= (l - |G| 2 ) 2vr 6(u + u') e (w) (18) 

This result does not depend on the specific amplification 
process. For a gain larger than unity, the added noise 
can be represented by a free field 6 m with the usual com- 
mutation relation (la) 



B in [u] 



i b" 



\J\G\ 2 - 1 (b ia [u]) 1 



(19) 



The presence of the conjugation is characteristic of am- 
plification processes and is encountered as soon as gain 
is present |2l]]2^] . 

We may use this example to describe the noise analysis 
in a measurement process. Let us consider that the field 
a ln carries a signal A superimposed with fluctuations c ln 



(20) 



The input noise power corresponds to the fluctua- 
tions a™ 



J AA 



(21) 



The measurement corresponds to the output a out of the 
amplifier 



= GA + Gc" 



\G\ 2 - 1 6 int 



(22) 



To analyze the noise of this amplified signal, we define 
an estimator A by normalizing the output field a out of the 
amplifier so that it is the sum of A and an extra noise 
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The added noise i s then described by a spectrum 

1 



(23) 
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(24) 



In the limit of large gain G and for thermal fluctuations, 
it corresponds to 



(25) 



When the temperatures are equal, this corresponds to a 
loss of 3dB in the signal to noise ratio. This effect has 
been observed since the beginning of radiowave commu- 
nications. It also sets a limit in the number of repeaters 
in optical fiber communications 23 241]. 



Most practical applications of amplifiers in measure- 
ments involve ideal operational amplifiers operating in 
the limits of infinite gain, infinite input impedance and 
null output impedance. In order to deal with the patholo- 
gies that could arise in such a system, we consider that 
it operates with a feedback loo p w hich fixes its effective 
gain and effective impedances p5[. 




a b 

FIG. 3. Representation of the ideal operational amplifier 
working in the limit of infinite gain with a reactive feedback 
Zf. (a) The noise sources are described as a current generator 
/ and a voltage generator U. (b) The amplifier is represented 
with a left (input) port I and a right (output) port r and the 
noise sources are modeled as input fields in the two noise lines 
a and a' 

We first analyze the amplifier as depicted on figure ||a 
where the noise sources are represented as a current gen- 
erator / and a voltage generator U. By coupling two 
coaxial lines denoted I and r respectively on the left port 
and the right port of the amplifier, one realizes a mea- 
surement model. The left line comes from a monitored 
electrical system so that the inward field l ln plays the 
role of the signal to be measured. Meanwhile, the right 
line goes to an electrical meter so that the outward field 
r out plays the role of the meter readout. In connection 
with the discussions of Quantum Non Demolition mea- 
surements pt|p7| , l out appears as the back-action field 
sent back to the monitored system and r ln represents the 
fluctuations coming from the readout line. A reactive 
impedance Zf acts as feedback for the amplifier. 

We now present the electrical equations associated 
with this measurement device. We first write the char- 
acteristic relations between the voltages and currents 



U 
I 



Ul = U r 
Il+If 



Z f I 



f 1 ! 



(26) 



Here, U p and I p are the voltage and current at the port 
p, i.e. at the end of the line p — I or r, while U and / 
are the voltage and current noise generators associated 
with the operational amplifier itself (see Fig.l). Zf is the 
impedance feedback. All equations are implicitly writ- 
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ten in the frequency representation and the impedances 
are functions of frequency. Equations (|26| ) take a simple 
form because of the limits of infinite gain, infinite input 
impedance and null output impedance assumed for the 
ideal operational amplifier. We also suppose that the 
fields incoming through the various ports are uncorre- 
cted with each other as well as with amplifier noises. 

As already emphasized, the output fields p out obey 
the commutation relations (^) of free fields. To make 
this property explicit, we use the characteristic equations 
(^,^2| ) associated with the amplifier and the lines to 
rewrite the output fields Z out and r out in terms of input 
fields l ln , r ln and of amplifier noise sources U and I 



-r + 



_ h\u\Ri 



u 



(27) 



We then deduce from ( |27| ) that the voltage and current 
fluctuations U and / obey the following commutation 
relations 



[t/M,E/[a/]] =[JM,I[a/]]=0 
[U [uj] , I [uj']\ = 2tt fiuj S (to + lo') 



(28) 



Hence, voltage and current fluctuations verify Heisenberg 
inequalities which determine the ultimate performance 
of the ideal operational amplifier used as a measurement 
device pEf . 

To push this analysis further it is worth introducing 
new quantities a ln and a' m as linear combinations of the 
noises U and / depending on a factor R having the di- 
mension of an impedance 



U M = y/2h\w\R (a in [u] - a' in [- 
I[uj] ~- 



2h\ 



R 



(a m [uj} + a nn [-u]) 



(29) 



For an arbitrary value of R, the quantities a ln and a'™ 
satisfy the free field commutation relations. In other 
words, the voltage and current noises associated with the 
amplifier may be replaced by the coupling to 2 further 
lines a and a' and the presence of amplification requires 
a conjugation of fluctuations coming in one of these two 
lines. This representation of the amplifier as a quantum 
network is depicted on figure |^b. 

We may then fix the parameter R to a value R a chosen 
so that the fluctuations a ln and a' m are uncorrelated. 
This specific value is determined by the ratio between 
voltage and current noise spectra 



R a = 



cru 



(30) 



The 2 noise spectra ajju and cru are defined as symmet- 
ric correlation functions. The fields a ln and a' m are thus 
described by temperatures T a and T Q / . We have assumed 
that these fluctuations are the same for all field quadra- 
tures, i.e. that the amplifier noises are phase-insensitive. 
Although this assumption is not mandatory for the forth- 
coming analysis, we also consider for simplicity that the 
specific impedance R a is constant over the spectral do- 
main of interest. 



IV. THE COLD DAMPED ACCELEROMETER 

We come to the discussion of the ultimate performance 
of the cold damped capacitive accelerometer designed for 
fundamental physics experiments in space p6[. 
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FIG. 4. Scheme of the capacitive sensor. The proof mass 
is placed between two electrodes formed by the inner walls 
of the accelerometer cage. The position dependent capaci- 
tances are polarized by an AC sinewave source which induces 
a mean current at frequency cut in the symmetrical mode. The 
mass displacement is read as the current induced in the an- 
tisymmetric mode. An additional capacitance is inserted to 
make the antisymmetric mode resonant with io t . The signal 
is detected after an ideal operational amplifier with capaci- 
tive feedback followed by a synchronous demodulation. The 
impedance of the detection line plays the role of a further 
resistance R r . The detected signal then feds the servo loop 
used to keep the mass centered with respect to the cage. 

The central element of the capacitive accelerometer is 
a parallelepipedic proof mass placed inside a box. The 
walls of these box are electrodes distant from the mass 
off a hundred micrometers. The proof mass is kept at the 
center of the cage by an electrostatic suspension. Since a 
three dimensional electrostatic suspension is instable, it 
is necessary to use an active suspension. 

In the cage reference frame, an acceleration is trans- 
formed in an inertial force acting on the proof mass. The 
force necessary to compensate this inertial force is mea- 
sured. In fact, as in most ultrasensitive measurements, 
the detected signal is the error signal used to compensate 
the effect of the measured phenomenon. 

The essential elements of the accelerometer are pre- 
sented in figure |^. The proof mass and the cage form 
two condensators. Any mass motion unbalances the dif- 
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ferential detection bridge and provides the error signal. 
In order to avoid low frequency electrical noise, the elec- 
trical circuit is polarized with an AC voltage with a fre- 
quency of a hundred kilohertz. After demodulation, this 
signal is used for detection and as an error signal for a 
servo control loop which allows to keep the mass centered 
in its cage. 

Furthermore, the derivative of this signal provides a 
force proportional to the mass velocity and simulates a 
friction force. This active friction is called cold damping 
since it may be noiseless. More precisely, the effective 
temperature of the fluctuations of this active friction is 
much lower than the physical temperature of the device. 

The detection is performed with the output detection 
signal r° ut . It is a linear combination of the external 
force F ext and of input fields in the various noise lines. 
We normalize this expression so that the coefficient of 
proportionality appearing in front of the external force 
F ext is reduced to unity. With this normalization, we 
obtain a force estimator F ext which is just the sum of 
the true force F ext to be measured and of an equivalent 
input force noise. In the absence of feedback, the force 
estimator reads flltj: 



Ma"" 



(31) 



where a 111 denote the various input fields corresponding 
to the active and passive elements in the accelerometer. 

When the feedback is active, the servo loop efficiently 
maintains the mass at its equilibrium position and the ve- 
locity is no longer affected by the external force F ext . The 
residual motion is interpreted as the difference between 
the real velocity of the mass and the velocity measured by 
the sensor. This means that the servo loop efficiently cor- 
rects the motion of the mass except for the sensing error. 
However the sensitivity to external force is still present 
in the correction signal. Quite remarkably, in the limit of 
an infinite loop gain and with the same approximations 
as above, the expression of the force estimator F ext is the 
same as in the free case Jl6| . 

The added noise spectrum £ff is obtained as 



of the electrical detection circuit. R a and 9 a are the 
characteristic impedance and temperature of the ampli- 
fier. 

In these conditions, the added noise spectrum is dom- 
inated by the mechanical Langevin forces 

= 1.1 x 1(T 25 (kg m s~ 2 ) 2 /Hz (33) 
This corresponds to a sensitivity in acceleration 



FF 



M 



= 1.2 x 



1(T 12 m s -2 



/VHz 



(34) 



Taking into account the integration time of the experi- 
ment, this leads to the expected instrument performance 
corresponding to a test accuracy of 10~ 15 . 

In the present state-of-the-art instrument, the sensitiv- 
ity is thus limited by the residual mechanical Langevin 
forces. The latter are due to the damping processes in the 
gold wire used to keep the proof mass at zero voltage || . 
With such a configuration, the detection noise is not a 
limiting factor. This is a remarkable result in a situation 
where the effective damping induced through the servo 
loop is much more efficient than the passive mechanical 
damping. This confirms the considerable interest of the 
cold damping technique for high sensitivity measurement 
devices. 

Future fundamental physics missions in space will re- 
quire even better sensitivities. To this aim, the wire will 
be removed and the charge of the test mass will be con- 
trolled by other means, for example UV photoemission. 
The mechanical Langevin noise will no longer be a limi- 
tation so that the analysis of the ultimate detection noise 
will become crucial for the optimization of the instrument 
performance. This also means that the electromechanical 
design configuration will have to be reoptimized taking 
into account the various noise sources associated with 
detection [EH. 



Tjff — 



E 
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We have evaluated the whole noise spectrum Y,p f for the 
specific case of the instrument proposed for the /zSCOPE 
space mission devoted to the test of the equivalence prin- 
ciple. Some of the main parameters of this system are 
listed below 



M = 0.27 kg 
i 5 x 1CT 4 Hz 



o 

2^ 

R a = 0.15 x io 6 n 



H m = 1.3 x 10" 

^ ~ 10 5 Hz 
2tt 

e Q = 1.5 k 



kg s" 1 



(32) 

M is the mass of the proof mass, H m is the residual me- 
chanical damping force, ^ is the frequency of the mea- 
sured mechanical motion, |^ is the operating frequency 
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